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ABSTRACT 

Context. Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 pm (called ’’coreshine”) seems 
to indicate the presence of grains with sizes above 0.5 pm. 

Aims. We aim to analyze 3.6 and 4.5 pm coreshine from molecular cloud cores to probe the largest grains in the size distribution. 
Methods. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 pm bands for a sample of 10 low-mass 
cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface 
brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. 
We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. 

Results. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes 
around 0.9, 0.5, 0.65, 1.5, 0.6, and ^1.5 pm, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the 
individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic 
local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho 
Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. 
Conclusions. Considering surface brightness ratios in the 3.6 and 4.5 pm bands across a molecular cloud core is an effective method 
of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth. 
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1. Introduction 

Dust grains are present in most cosmic objects, and they add a 
rich spect rum of investigation methods to astrophysical research 
(see, e.g.. lDrainell2003l:lHenningll20 1 Ql) . Through the absorption, 
emission, and scattering of radiation they allow us to trace the 
local density and temperat ure structure of objec t s (for molec¬ 
ular cloud cores see, e.g.. iNielbock et al.l I2Q12I: iLippok et al.l 
I2Q13I: lUaunhardt et~^l2Q13l) . and to reveal the action of young 
stellar object s (YSOs), which are enshrouded at short wave¬ 
lengths (e.g.. iHenning et al1ll99Ql) . In the interstellar medium 
(IS M), dense regions are predominantly cooled by dust radia¬ 
tion (lGoldsmithlf200ll) . Chemical reactions o n their surface con- 
tribute to the chemical processes in the gas (iHerbst et al.ll2005l) 
with the available surface depending on the grain size distribu¬ 
tion and the fractal degree. Dust grains are sens itive to magneti c 
forces because of their charge and composition (lLazarianll2007l) . 
In the planet formation process, small dust grains are the seeds 
for forming larger bodi es in the accretion disks around YSOs 
(ISteinacker et al.|[2013ah . 

Important questions therefore include the origin of the 
grains, their composition, and size distribution. Dust grains 


are k nown to be produced and destroyed in supern ovae events 
(e.g., IWesson et all I201QI: lindebetouw et all l2014l) and in the 
atmospheres of stars in their late evolutionary phases (e.g., 
[McDonald et al.|[200^ . Their distribution into the ISM and their 
complex processing includ ing destruction and g r owth is the sub¬ 
ject of ongoing research (Uones & l2Qlll: lAndersen et al.l 
[ 2011 !: IZhukovska & Hennin gl |20 1 4l) . Grain growth processes 
have been discussed mostly in the framework of accretion disks 
and planet formation. Cold and dense prestellar molecular cloud 
cores may also exist long enough t o allow the collisional growth 
of grains (e.g.. (IOssenkopj|l993l:IOssenkopf & Henning 1994 


01 grams (e.g.. (IUssenkoptlli99JI:IUssenkopi & LLenmng 1199^: 

IWeidenschilling & Ruzmaikinalll994l0rmel et al.ll2Q09L 201 ll) . 


There is evidence t hat grains are ice-coated in the dense ISM 
(IWhittet et al.ll 19831) . supporting the growth process by efficient 
sticking. However, it remains to be investigated to what extent 
gas turbulence can provide the relative velocities that are needed 
to coagulate grains as assumed in these molecular cloud core 
models. 


Co ncerning modeled grain size distributions, iMathis et~^ 
(I 1977 I) were the first to propose a distribution for the diffuse 
ISM. They suggested a power-law distribution with a slope of 
-3.5 and with a size limit of 0.25 pm for non-graphite grains 
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(abbreviated as ”MRN distribution”) and emphasized that only 
weak constraints could be derived for the larger grains in the dis¬ 
tribution. In the following we denote size distributions as ”MRN- 
type” when they have a power-law index of -3.5. When we refer 
to grain sizes, we will give the grain radius value. Since then 
a large number of observations and models considering thermal 
emission and extinction of grains have suggested both varying 
size distributions and gra ins with radii beyond 0.25 jjm 'm the 


denser parts of the ISM (ISteonik et al.ll2003l: iKiss et al. 

2006 

Ridderstad et al .11 20061: ISchnee et al.ll2008l:IChaDman et al. 

2009 

Chapman & Mundvl 

20091: iMcClurel 120091: iParadis et all 

2010 

Veneziani et al.ll2010 

: Ysard et al.ll2013l). 


Scattering by dust grains is another process that is sensi¬ 
tive to the grain size. iLehtinen & Mattila (Il996h modeled the 
1.25, 1.6, and 2.2 yum band (J, H, and Kn) scattered light im¬ 
ages of the Thumbnail nebula and concluded that the limiting 
grain size exceeds 0.25 //m. The strong decline of scattering ef¬ 
ficiency of ISM grains beyond the Ks band caused people to 
assume that the cores become dark in scattered light in the mid- 
infrared (MIR). It would take even larger grains of micron-size 
scale to elevate scattered light intensities back to the surface 
brightness (SFB) measurable by MIR detectors. Such radiation 
was indeed found insid e the cloud LI83 and named coreshine 
(ISteinacker et al.ll201Ql). This t erm di stinguished it from cloud- 
shine, as [Foster & Goodman! (l2006l) named the near-infrared 
(NIR) scattered light from the core parts that are not shielded 
from NIR radiation. Beneficial for the detection of this radiation 
was the fact that extinction has a minimum in the wavelength 
range 3-5 //m, allowing interstellar radiation to penetrate deeper 
int o the core. _ 

ISteinacker et all (120101) present the most complex spatial 
structure modeling that has been performed so far for a molecu¬ 
lar cloud core. Based on SpitzerjlKAC data at 3.6, 4.5, and 8 yum, 
the 3D radiative transfer (RT) modeling of LI 83 required the as¬ 
sumption of population of grains with sizes up to 1 yum. They 
argue that the excess of SFB in the 4.5 yum band could not be 
fit by emission of transiently-heated grains due to the absence 
of a strong feature in this band. Despite the strong central ex¬ 
tinction (Ay above 150) and high Galactic latitude {b = 36.8°), 
the physical conditions in the LI83 clump are not intrinsically 
different from those seen in nearby cores. Correspondingly, the 
investigation data of Spitzer PI and Legacy survey data of a sam¬ 
ple of 110 cores led to th e conclusion that ab out half of all cores 
show hints of coreshine (jPagani et al.ll2010l3) . It was shown that 
scattered light from cores is observable in a variety of star forma¬ 
tion stages, from prestellar cores to disks around forming stellar 
objects, as well as for a wide range of morphologies, i.e. sin- 
gle or binary sy stems and in filamentary structures (see also 
IStutz et al] l2009l for a small subsample of cores with 3.6 pm 
emission). In the same study, the RT analysis of model cores re¬ 
vealed that without background, the scattered light morphology 
does not vary strongly with the position in the Galaxy aside from 
weak enhancement towards the Galactic center (GC). For more 
massive cores with 10 M©, the morphology showed a stronger 
crescent of enhanced emission due to the combined action of 
extinction and multiple scattering. 

Since then further studies have investigated the occurrence 
and conditions of coreshine detection, and this work is also de¬ 
voted to scattered MIR r adiation in niolecula r clouds and the 
grain properties causing it. iPagani et~n] (l2Q12l) find a drastically 
reduced occurrence of scattered light (3 detections plus 3 uncer¬ 
tain cases out of 24 objects) in the Gum/Vela region. It has been 
suggested that a connection to the action of a supernova rem¬ 
nant blast wave affects the grain size population in this region. 


In general the analysis of coreshine in regions with an enhanced 
interstellar radiation field (ISRF) like the Vela region is impeded 
owing to the confusion with emission of stochastically-heated 
grains like PAHs (polycyclic aromatic hydrocarbons). 

Observations in th e bands at 2.2, 3.6, an d 8 pm of the core 
L260 were modeled bv lAndersen et al.l(l2Q13l) . Their simple core 
model was able to reproduce the observed SFB profiles that as¬ 
sume an MRN-type size distribution with grains up to about 1 
yum in size and a locally enhanced incident radiation field (en¬ 
hancement factor 1.7). 

ISteinacker et all (l2Q14bl) have analyzed the coreshine seen 
from the low-density core L1506C and find that it requires the 
presence of grains with sizes exceeding the MRN distribution. 
Using grain growth models, they argued that L1506C must have 
passed through a period of higher density and stronger turbu¬ 
lence to grow the grains in situ. 

The conditions of detecting scattered light from low- 
mass molecular cores a t 3.6 yum have been re-investigated in 
ISteinacker et al.l (l2014al) including background effects. The lim¬ 
its derived from the RT equation indicate that extinction by the 
core prohibits detection in bright parts of the Galactic plane and 
especially near the GC. They show that the scattering phase 
function favors the detection of scattered light above and be¬ 
low the GC and, to some extent, near the Galactic anti-center. It 
was discussed that an enhanced radiation field may give rise to a 
coreshine signal even in the presence of a strong background. 

iLefevre et all (l2014l) have investigated the key parameters 
for reproducing the general trend of SFBes and intensity ratios of 
both coreshine and near-infrared observations. Based on a care¬ 
ful determination of the background field, they find that for a 
sample of 72 sources, starless cores show 3.6yum/4.5 pm SFB ra¬ 
tios above 2, while cores with embedded sources can have lower 
values. To constrain the dust properties, they used a rich grid 
of size-averaged dust models based on an extrapolation of stan¬ 
dard spherical grains able to fit the observations in the diffuse 
medium. Moreover, the effects of fiuffiness, ices, and a handful 
of classical grain size distributions were also tested. As spatial 
density distribution, an inclined ellipsoid with a Plummer-like 
profile with two masses was assumed. Their density structure 
is well-suited to centrally-condensed cores like LI544, which 
show a central depression in their SFB pattern, while the vast 
majority of cores with coreshine have no depression. They find 
that normal interstellar radiation field conditions are sufficient to 
explain coreshine with suitable grain models at all wavelengths 
for starless cores. According to their multi wavelength approach, 
the standard interstellar grains are not able to reproduce obser¬ 
vations and only a few grain types meet the criteria set by the 
data. 

This paper investigates coreshine in a sample of ten cores. 
Coreshine has been detected for all cores except for Rho Oph 
9 where the emission could also be explained by expected PAH 
emission. We use deep warm Spitzer IRAC 3.6 and 4.5 yum band 
observations (Cycle 9 Coreshine Follow-up: Program ID 90109, 
PI R. Paladini, Paladini et al., in prep.)Q . The observations were 
motivated by the fact that silicate grains of a size around 1 pm are 
expected to have their largest scattering-to-absorption efficiency 


^ ’’Cold” and ’’Warm” denote, respectively, the first part of the Spitzer 
mission - which lasted from August 2003 to May 2009 - during which 
all the instruments, including IRAC, were cooled to cryogenic temper¬ 
atures, and the second part of the mission - from May 2009 to present 
- which started at the end of cryogen and in which only the channels 1 
and 2 of the IRAC instrument are operating. 
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in the wavelength range 3-5 yum and that the largest grains should 
have their strongest impact in the 4.5 yum band. 

In Sect. 2 we briefly summarize the data and their acquisi¬ 
tion and processing. The SFB ratio model is derived in Sect. 3, 
outlining the different RT transfer methods and optical depth ef¬ 
fects. The modeling is described in Sect. 4, and the maximum 
grain sizes are analyzed for each core assuming an MRN-type 
size distribution. In Sect. 5, we discuss the results and consider 
the influence of local isotropic radiation fields, and summarize 
our findings in Sect. 6. The determination of the off-core SFB 
is described in Appendix A, and in Appendix B we describe the 
opacities of the dust model. 

2. Data 

The 3.6 and 4.5 yum IRAC data used for this work were taken 
during cycle 9 of the Spitzer warm mission. These observa¬ 
tions targeted a set of ten Galactic cores partially selected from 
the ’’Hunting for Coreshine” program (Program ID 80053, PI 
Paladini). In earlier data, nine cores show a clear indication 
of coreshine emission at 3.6 pm. Extinction at 8 pm can be 
found for all cores in the cold Spitzer data except for PLCKECC 
G303.09-16.04 (we use the abbreviation ecc806 in this paper) 
and L1506C. It is a core from the Planck Early Release Cold 
Core Catalogue (ECC), which is part of the Planck Early Release 
Compa ct Source Catalog (ERCSC) (iPlanck Collaboration et all 
We selected a varied sample of cores for this study, 
some with particularly strong coreshine (e.g., L260, ecc806, 
L1262, L1512), some also previously well-studied for coreshine 
emission (L260, L1506C) or otherwise (L1544, L1498). Some 
of these and other cores have particular features of interest: bi¬ 
nary core L1262, L1517A having similar cores nearby, L1512 
with a particularly simple shape, LI544 with a central coreshine 
depression, L1506C and L1498 with unusually low densities, 
L1439 near more evolved young stellar objects, and Rho Oph9 
with an expected strong PAH contribution. The original survey 
did not reach a comparable signal-to-noise (S/N) at both 3.6 and 
4.5 yum, with the consequence that coreshine emission at 4.5 yum 
was typically at the level or even below the instrumental noise. 
This situation was rectified in Cycle 9. These new observations 
were carried out in full array mode with 100 s integration. A 
detailed description of the source selection, observing strategy, 
and data processing will be provided in a forthcoming paper 
(Paladini et al., in prep.). In this work we have used the IRAC 
Level 2 mosaics (or pbcd, post basic calibration data), generated 
by the Spitzer Science Center (SSC) pipeline and downloaded 
from the Spitzer Heritage Archive (SHA). These images are af¬ 
fected by column-pulldown, an instrumental effect that consists 
in a depression of the intensity registered by pixels along a given 
array column, and which is caused by a bright source, i.e. a star 
or a cosmic ray. The IRAC images are also not zody-subtracted, 
meaning that the contribution from zodiacal light is not removed 
by individual frames or from the final mosaics at any stage of the 
automatic processing. 

Basic core properties gathered from the literature are given in 
Table [T] The sources are ordered by decreasing coreshine peak 
SEB (see also Table O. The estimates for the central column 
density have error bars (factors of a few) partly because the val¬ 
ues are derived using a specific dust model. The ranges given in 
the Table are either specified in the cited work or based on as¬ 
suming uncertainties by a factor of 2. Values without citation are 
based on simple estimates and just serve as an indication. The 
given masses and sizes (mean core diameter) are derived using 
different methods and also depend on the definitions of the core 


boundary. Moreover, it usually takes more than one parameter to 
characterize the complex core structure. In some cases we have 
estimated the size from the Spitzer map. To overcome the varia¬ 
tion in the optical depth due to the uncertainties, we will rely in 
the modeling on a maximum optical depth at 3.6 pm (2.2 pm in 
the case of L260) that is given in the last column. Its choice and 
role in the modeling is described in Sect. 13.71 Eor modeling in¬ 
dividual cores with more sophisticated density models and using 
multiwavelength 3D RT, it will be worth re-evaluating the mass 
estimates based on the original data and using a more sophisti¬ 
cated structure model with more than just a few size parameters. 
The meaning of Tmax is explained in Sect. 3.7. 

To give an impression of the overall SEB levels, the source 
geometry, sizes, and the influence of neighboring stellar sources, 
we show in Eig.[T]the 3.6 yum off-subtracted SEB images of all 
sources scaled to the same distance. Eor the example of L260, 
Appendix A describes how we have determined for all cores 
this ”off-core surface brightness” that is subtracted from the 
core SEB. The linear scale in the plane-of-sky (PoSky) is in¬ 
dicated by the 10 kau bar in the lower lefthand corner. To in¬ 
dicate the approximate location of the central part of the core, 
we have marked the minimum of the corresponding 8 pm maps 
(for L1498, L1262, L260, L1439, L1544, L1517A, and L1512) 
where the cores are seen in extinction or the maximum of ther¬ 
mal emission maps in the EIR/mm (ecc806 and L1506C). Details 
are given in Sect. 4. 


3. Modeling approach 

Coreshine modeling involves RT calculations in a complex 3D 
core geometry. The scattering integral in the RT equation mixes 
radiation from all directions in all points. Moreover, to describe 
a 3D density structure, about ten to hundreds of free parame¬ 
ters are needed depending on the complexity of the structure. 
This makes coreshine modeling with an anisotropic illumination 
of the dust grains numerically difficult. A scan of the entire pa¬ 
rameter space of size distributions, chemical compositions, grain 
shapes, and spatial density distributions is therefore too demand¬ 
ing for the current computer power. In the following sections, we 
describe our basic assumptions for reducing the parameter space 
and the complexity of the problem, the RT approximations, and 
the general modeling approach. 

3.1. Basic assumptions 

To reduce the complexity of the problem we make three basic 
assumptions. 

Assumption 1: The grain size distribution s{a) da is normalized 
to have the same gas-to-dust mass ratio across the core, and it 
has no spatial change in shape or size limits. 

This assumption is supported by the shallow gradients of the 
SEB pattern observed across the cores with coreshine. In turn, 
ISteinacker et all (l201Qh used a grain model with grain size in¬ 
creasing with gas density, and the SEB model showed stronger 
gradients than observed. Nevertheless, there are theoretical ar¬ 
guments for why the opacities might change across the core. Eor 
example, grain coagulation and the thickness of the ice mantles 
are expected to increase toward the center of the cores, depend¬ 
ing, however, on how quic kly turbulent mixing c an smooth these 
gradients again (see also I Andersen et aHl2013l) . Relaxation of 
this approximation will be the subject of further publications. 
Assumption 2: The size distribution is a power law, and the sizes 
range from ai io 02 . 
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Fig. 1. Off-subtracted surface brightness images at 3.6 yum for the sources of the sample. The image sizes have been adjusted to 
show the sources at the same distance in order to visualize the source sizes. Most cores have a complex shape, and the diffuse 
signals differ with respect to the peak brightness and point source contamination. In the case of Rho Oph 9, the image has not been 
off-subtracted and the bright emission in the upper left corner is caused by transiently-heated particles and not by coreshine. The 
crosses indicate the approximate location of the core center based on the corresponding 8 //m extinction maps or thermal emission 
maps in the FIR/mm. 


We have varied the power-law index and discuss its impact 
on the findings in Sect. 14.11 While there are studies suggesting 
that the spectral index of the d istribution may be different fm m 
-3.5 as observed in the ISM (IWeingartner & Drainell2QQll) or 
deviate from powe r laws (IWeidenschilling & Ruzmaikinall 19941: 
lOrmel etd] l2009l) . we use this distribution as a typical case 
where the total grain surface area is dominated by the small 
grains while the mass predominantly rests in the large grains. 
It needs to be verified whether currently available data allow 
constraining the detailed size-distribution beyond the power-law 
picture. 


Assumption 3: The optical properties of the grains are derived 
from Mie theory using spherical ice-coated silicate and carbona¬ 
ceous grains. 


Numerical calculations by (lOrmel et al.ll2Q09l) that are based 
on aggregate collisions show that grains in cores go through a 
phase of compaction so that the spherical approximation might 
be better than initially expected from a growth picture that 
involves fluffy aggregates with large geometrical Ailing fac¬ 
tors. Nevertheless, future modeling with a more complex dust 
model should inclu d e non-spherical grain s (as discussed, e.g., in 
iLefevre et al.ll2014l) . IWhittet ef^] (119831) showed that grains in 
the Taurus cloud contain ice , and later stud i es hav e confirmed 
this finding for other clouds. lAndersen et al.l (l2014l) And for the 
Lupus IV molecular cloud complex that the limits for the oc¬ 
currence of ice mantles and coreshine are similar and possibly 
related. We use a mass ratio of 1:4 for the carbonaceous and 


4 





























Steinacker et aL: Grain size limits derived from coreshine 


Table 1. Core sample and their basic properties'* 


Source 

Gal. long. 

n 

Gal. lat. 

n 

Distance 

[pc] 

N. 

[10^'^ m-2] 

Mass 

[Mol 

Size 

[kau] 

C max 

at 3.6 pm 

L260 

008.65 

+22.\7 

160^ 

1-4 


36" 

1 at Ks"* 

ecc806^ 

303.06 

-16.04 

160^ 

0.5-2^ 

2-8« 

25" 

1 

L1262 (CB244) 

117.12 

-h 12.40 

200^ 

3-21P 

3-T 

15*^ 

4.5 

L1517A 

172.50 

-07.98 

140" 

1-r 

L 

20" 

1 

L1512 (CB27) 

171.87 

-05.24 

140" 

1 - 9 P 

2.4^ 

15" 

1.5 

L1544 

177.98 

-09.72 

140" 

6-26^ 

2.8^ 

20^ 

6 

L1506C 

171.15 

-17.57 

140" 

1-4 

4h 

25" 

0.5 

L1439 (CB26) 

156.06 

+06.00 

140^ 

0.1-6P 

1.6^ 

15^ 

0.5 

L1498 

170.14 

-19.11 

140" 

1.6-6.4^ 

0.51-0.83*' 

30^' 

0.5 


References. IVisser et al.l ( 200l). ^Caselli et al. 120021). ^read from th e Spitzer image using t h e distance. IStutz et alJ d2 01Ql). iKenyon et al.l(ll994l) . 
iHacar & Tafallal (1201 ll) . ^ ~utz et alJ (120091) . ^ Pagani et all (l2010ah . IShirlev et~Sl j2005h . ^^fonaetani^OO^r^aseffije^hTloOS^Tiofficml 
name is PLCKECC G303.09-16.04, "*for L260, we used a constraint in the 2.2 //m band Ks, "^derived from iHacar & Tafallal (1201 Ih. dLippok et al 
2013h . ^based on Herschel/SPIRE archive data of the Gould Belt Survey (PI: Andre) '^Rho Oph 9 is not included. iPlanck Collaboration et al 


silicate component and assume the same size distribution. The 
details of the dust model are described in Appendix |B1 


3.2. Optical depth 


The essential quantity for characterizing the impact of the dust 
on the transport of radiation is the size-integrated optical depth 
between two points on a line x = Xq £n with the starting point 
jco and the direction n. Based on the approximations listed in 
Sect. l3.1l it has the form 


(r), (T, A, ^1,^2) 


i\ a2 



da s(a) cr(A, a) 


0 ai 


( 1 ) 


where t is the spatial distribution of the H 2 number density, and 
A the wavelength. Integrating t over space gives the total number 
of H 2 and, when multiplied with the H 2 mass, the assumed total 
core mass. Since (r)^ contains the line-of-sight (LoS) integral 
of the gas density, it can be related to the H 2 number column 
density N between the two points. The cross section cr refers 
to absorption, scattering, or extinction, which defines the optical 
depth (Tabs)s^ (rsca)s^ and respectively. 

An important feature of the optical depth is that it depends 
not only on the core properties but also on the dust properties, 
so that for the same core, for example, two dust models can give 
very different central optical depths. This coupling of the core 
density structure and the dust properties with the optical depth is 
not easy to disentangle: basic properties like total core mass, its 
size, or visual extinction are often derived by interpreting obser¬ 
vational data with a specific dust model and will change when 
the modeling requires a change in the dust properties. 

In this paper, we use three approximations of the general 3D 
RT equation for this application based on the assumed maximum 
optical depth and give them simple names with the precise defi¬ 
nitions provided in this section. 


3.3. Tull RT 

For cores reaching optical depths for extinction beyond I, shad¬ 
owing effects and multiple scattering become important in the 
core. This is the most complex of the three modeling variants, 
and the stationary 3D RT equation in this case has the form 

dl{x, n, A) 

-TT-= - {crext)s Cii,a2) t{x) I(x, n,A)-\-S (x, n, A)-\- 


-\-t(x) J" dO! {(Ts, 

An 


, p), (T, fi, Q', < 21 , < 22 ) /(x, n\A). 


( 2 ) 


The equation describes the change in the intensity I given at the 
location x in the direction n along the path dl. The source term 
S contains the incident interstellar radiation field or nearby radi¬ 
ation sources. We have abbreviated the integrals over the grain 
sizes by 





da (Texti^, u) s{a) 


( 3 ) 


and 


{CTsca P)s 



da CTscai'^^ cl) p(A, a, n) s(a). 


(4) 


The probability of scattering radiation from the solid angle JQ' 
into the considered direction n is given by the phase function 
p. This non-local integro-differential equation requires applying 
advanced numerical schemes like Monte Carlo or ray-tracing in 
order to solve for the inte nsity emerging in the d irection of the 
observer (for a review see lSteinacker et al.]l2QI3bl) . 

An essential difficulty in analyzing radiation received from 
a core with optical depths beyond I is that the structural in¬ 
formation is mixed in at all LoS both by extinction effects and 
by multiple scattering. Since most cores show complex spatial 
structures, the error from assuming a wrong density model will 
be propagated with every ray or photon propagation into the de¬ 
rived SFB. It is important to note that this is not a problem of the 
available solvers. It is a problem of the modeling, which tries to 
avoid using too many free parameters that enter with the asym¬ 
metric core structure. Shadowing effects also amplify the impact 
of anisotropic illumination. The RT modeling of cor eshine from 
the core LI 83 performed in I Steinacker et al.l (l20IQl) . for exam¬ 
ple, required the use of 100 3D Gaussian density clumps. An 
illustration of the RT variants is given in Fig. [2] with the full RT 
case being described in the righthand panel. 


3.4. ’’Single-scattering RT” 

When the optical depth for scattering is below I across the core, 
the impact of photons that scatter more than once is reduced. In 
this case we can use the single-scattering approximation. Instead 
of following rays or photons in arbitrary directions, we can sep¬ 
arate contributions along the LoS (x-direction) in each PoSky 
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Fig. 2. Illustration of the different RT methods used in the paper (the abbreviation ’’Thin RT” stands for ”RT in the optically thin 
limit”). 


position (y,z). Constructing the SFB from the background in¬ 
tensity It,g that undergoes the extinction and radiation that is 
scattered into the LoS from all directions, we find (see, e.g., 
ISteinacker et al.ll^l4ah 


I(A,y,z,ai,a2) = e 

+00 aj 


= + Ifg(A)+ 

J' dx t(x) J' da s{a) CTscaid, a) {IP}q (a, x, y, z) 


(5) 


with the foreground SFB the direction-averaged product of 
incident radiation field, and phase function P now giving the 
probability of scattering radiation from the solid angle JQ into 
the LoS 


is a result of the combined action of opacity and spatial density 
effects. 

3.5. RT in the optically thin limit (’’Thin RT”) 

For some cores with low column densities and/or low-opacity 
grains, it is possible to treat RT in the single scattering and op¬ 
tically thin approximation even for grains with sizes of several 
pm in the wavelength range 3-5 pm. The limiting r values for 
the thin approximation can be related, for example, to the con- 
sidered order of the Tayl or expansion of the exponential function 
(ISteinacker et al.|[^14ah . 

The SFB for an expansion to first order in Text is 


{IP}q (a, X, y,Z) = e '^ext-cell-obsiAx,y,z) 

An 

Interstellar radiation on its way to the core position (x, y, z) 
from a certain solid angle Q undergoes extinction described by 
^ext\sky-ceih The radiation that is scattered by the dust at (x,y,z) 
and leaves towards the observer again undergoes extinction with 
the optical depth Text-ceii-obs- 

The extinction of radiation still carries 3D density informa¬ 
tion from the entire core into the intensity of radiation being scat¬ 
tered into the LoS, and with it all errors of the density approx¬ 
imation. But the directions of the radiation transport are known 
and can be analyzed in advance. 

The advantage of this approach is that it allows precalcu¬ 
lation of all number column densities that are encountered by 
photons arriving from all directions at the considered LoS and 
then propagate to the observer. This speeds up the calculations 
when several wavelengths are considered and allows for the use 
of higher resolution ISRF maps, for example. We have created 
a code version that m akes use of this advantage (used, e.g., in 
I Andersen et al.ll2Ql^ and apply it also in this paper. 

Owing to optical depth effects, the two bands at 3.6 and 4.5 
pm can probe different regions in the outer parts of cores with 
masses of a few M© when the opacities change from 3.6 to 4.5 
pm. For a standard core profile (fiat followed by power law with 
an index of -1.5), a mass of 5 M©, an outer radius of 10 kau, and 
an MRN-type distribution of silicate grains, among others, the 
locations of the r = 1 layers differ by 14 % of the core radius, 
and the probed column number densities differ by a factor of 1.8. 
For single-scattering RT like for full RT, the SFB in both bands 


I(A, y, z, au ^ 2 ) - 4^(4, y, z) - Ifg(A, y, z) ~ 
- {o-ext}s (T, au ^ 2 ) N(y, z) 4g(T, y, z) 


a2 

/ 


An 


+N(y,z) I da s{a) crsca{A, a) I dfl P(A,fl,a) Iin(A,fl). (7) 


The lefthand side of Eq. (0 and the radiation behind the core 
Ihg can be determined observationally. Also for the incident field 
Iin, an approximation can be derived from observed allsky maps. 
Then we can calculate the righthand side with a model for the 
density and dust opacities. Unlike the cases of full or single scat¬ 
tering RT, the spatial and dust property problem decouple for 
thin RT since both terms on the righthand side of Eq. (0 con¬ 
tain the spatial variation in the form of N for each pixel or LoS. 
Since for optically thin scattering, both bands see the same col¬ 
umn of grains, the spatial density variation is separated from the 
opacity variation. This makes the low-density cores best-suited 
to studying the optical properties of the grains as performed for 
translucent clouds and the diffuse ISM. 

3.6. Modeling surface brightness ratios instead of surface 
brightnesses 

Eor the modeling approach of the two images obtained at 3.6 and 
4.5 pm for each source, we considered three methods partially 
based on the ratio, R, of the off-subtracted images. 

i) Image modeling, where we assume a dust model and a 
density structure, calculate model images in the two bands, and 
compare model and real images at each wavelength. The model¬ 
ing process is to vary the density structure and opacity until the 
differences are minimized. 
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ii) '’Full RT” R and maximum surface brightness modeling, 
for which we build the ratio map from the observed images, se¬ 
lect a region where the optical depth is expected to be small, and 
determine the observed/? distribution. Then assuming a dust and 
density model, we calculate the model images in the two bands 
using ’’full RT”, build the ratio map, and find the theoretical R 
distribution in the selected region. From the images we find the 
maximum SFBs. The modeling process is to vary the opacity and 
the central column density until the differences are minimized. 

iii) ’’Thin RT” R and maximum surface brightness modeling, 
by building a ratio map R from the observed images in the two 
bands, selecting a region where the optical depth is expected to 
be small, and determining the observed R distribution. Then as¬ 
suming a dust model and a maximum column density, we can 
calculate the maximum SFB. The modeling process is to vary 
the opacity and the central column density value until we can 
reproduce the mean R and the maximum SFB. 

To avoid either very complex spatial modeling with many pa¬ 
rameters or the introduction of large errors from taking a density 
structure model that is too simple, we applied methods ii) and 
iii). To build the ratio map R = 4^(3.6 yum)//c5(4.5 jam) from 
the SFBs above the background caused by scattered light in the 
core Ics, we divided Eq. (0) for the two bands. For each pixel in 
the two images, the ratio on the lefthand side can be determined 
from the observed data 


Robsiy, z) 


r-Hy, z) - i^gHy, z) - if^iy, z)' 


( 8 ) 


Using the albedo 


w^(a, a\,a 2 ) = 






(9) 


and the mean product of albedo and phase function 


= J"da s(a) w^(a,ai,a 2 ) p^(Tl,a), (10) 


we can express the righthand side ratio for optically thin and 
single scattering as 


Rfheo^PX') ^ 2 ) 


/^3.6\ (O) - Ilf 

\_e^An - 

An 


This equation contains only the dust properties (cross sections, 
phase functions, and size distributions) and the incident field 
(all-sky map //„ and background 4^ of the core). 

Practically speaking, although some cores in the sample can 
be modeled with optically thin RT, we used single-scattering RT 
since the numerical effort is only marginally greater. The thin 
RT case is considered here because of its properties for separat¬ 
ing dust opacity and structure, an advantage that is kept to some 
extent for higher optical depths. 

To investigate the errors from using the three approximate 
RT schemes, we performed RT calculations for model cores with 
the masses typical of the range of low-mass cores 0.75, 1.7, 
and 2.4 M© and chose the core scale so that the optical depth 
for extinction through the center is 0.3, 0.7, and 1, respectively. 
Figure [3] compares thin and single-scattering RT in the upper 
panel. The images show the SFB error of the core at 3.6 pm 
\Ithin-Isingie\IImax in percent with the maximum SFB Imax of I full- 


As expected, the maximum error for r < 0.3 is small (about 2%). 
At T < 0.7, the error is in the 5% range, for r < 1 it reaches 12 
%. Comparing thin and full RT, the error up to r < 0.7 is the 
same since the extinction effects and multiple scattering can be 
neglected. For r < 1, both effects start to modify the resulting 
SFB, and the error rises to 15%. 


3. 7 . The ’’limited optical depth” approach 

To decide which core allows for which RT scheme, we need to 
specify the optical depth for the cores in the sample. However, 
neither the mass and density distribution in the core nor the op¬ 
tical properties of the dust grains are well known (the targets of 
investigation in this paper). The ranges in Table [T] give an indi¬ 
cation of uncertainties estimated by authors and are often factors 
of a few. 

Coreshine and scattered light from cores in general have a 
feature that can aid us in constraining the optical depth. When 
the central optical depth increases beyond one, the SFB profiles 
start to flatten and then show a central depression. The exact 
value when the depression is clearly visible depends on the prop¬ 
erties of the core, the incoming field, and the background field, 
but in our experience, t=2 is a reasonable value to expect the 
occurrence of depressions in the images. Figure [T] contains two 
prominent cases (LI544 and Rho Oph 9) and a less pronounced 
one (LI262), but most of the cores have a pattern with a cen¬ 
tral maximum. In the modeling, we can use this observational 
evidence to consider only models with a central optical depth 
below a few whe n no central depression is visible at the consid¬ 
ered wavelength. I Andersen et al.l (l2Q13l) used this argument to 
exclude grains with sizes beyond 1.5 pm since with the avail¬ 
ability of NIR data, it is possible to find the wavelengths where 
depression sets in. 

We therefore chose to apply a method which we call ’’lim¬ 
ited optical depth” approach. It adjusts density and opacities in 
a way that they do n ot exceed a given r limit. Fixing the den¬ 
sity (as done, e.g., in iLefevre et al .11201^ and then considering 
various dust models, r may rise above a few and produce a cen¬ 
tral depression while the data show a maximum in SFB near the 
projected core center. Especially, models with steep density pro¬ 
files and small kink radii (below which the profiles are assumed 
to be fiat) accumulate most mass in the center leading to high 
central column densities. Furthermore, central column densities 
are usually derived from thermal emission measurements in the 
FIR/mm wavelength range based on a specific dust model, us¬ 
ing the column densities with another dust model thus might be 
incorrect. The limited optical depth allows us to avoid this spec¬ 
ification and the dust properties are consistently used only in the 
modeling of the coreshine. 

In Table [T] we summarize the final maximum optical depth 
at 3.6 pm (2.2 pm in the case of L260) that we have used for 
the core modeling. Cores without central depression have values 
below 2. For the cores with Tmax > 2, we chose values that scale 
approximately with the average column density value given in 
the literature. The modeling is performed within the ranges of 
the mass and optical depth given in Table [T] 

3.8. Properties of the theoretical surface brightness ratio R 

In this section, we discuss how the cores can appear in the two 
bands, which basic physical properties are responsible for the 
appearance, and what the resulting range of theoretical SFB ra¬ 
tios is. Since the largest grains in the size distribution scatter the 
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Fig. 3. Relative error in using the RT approximations. 


radiation more efficiently, our goal is to explore the maximum 
grain sizes when comparing the observed R distribution across 
the core with the theoretical R of a particular dust model. 

Besides scattering the incident radiation toward the observer, 
the grains are also responsible for the extinction of background 
radiation along the LoS. This background extinction introduces 
several complications. The scattered light has to overcome the 
extinction of the background radiation by the core to arrive at a 
SFB above the background and to appear in excess of emission. 
The four possibilities are to see the core in extinction in both 
bands, in just one of them, or in emission in both bands. Since 
for our core sample we have analyzed only parts of the cores 
with excess of emission in both bands, we concentrate on the 
last case. 

To keep the example simple, we consider the condition for 
excess of emission in the optically thin case. From (fTTI) . we find 

< JdQ (wV)/^) ( 12 ) 

47r 

for both bands. The most prominent effect on R is the pole 
causing a singularity when the 4.5 jum denominator vanishes. 
Physically, this happens when the scattered light SFB is equal to 
the extincted background SFB. Increasing the maximum grain 
size, the SFB in the two bands increases due to the higher scat¬ 
tering efficiency of larger grains. In the picture of Rayleigh scat¬ 
tering, the condition 27ra T is less fulfilled for larger grains. 
The 4.5 jum band SFB benefits stronger from the increased max¬ 
imum size due to this condition and R decreases. 

Figure |4] shows in Panel A the numerator and denom¬ 
inator of R for the example of L260. We have used ice- 
coated silicate grains and an MRN-type size distribution rang- 
ing from 0.01 am to a 2 . We follow the approach described in 
I Andersen et al.1 (l2013l) to calculate the incoming fields from 


zodiacal-subtracted DIRB^ maps. The calculation is done in 
single-scattering approximation. 

The left yellow regions indicate where the core is seen in 
extinction in both bands or just the 4.5 jum band. The observed 
coreshine SFB defines the limit of this region when becoming 
zero. The thick dashed lines give the maximum observed core¬ 
shine SFB, and we also show the uncertainty in the value due 
to the choice of background (see Appendix O by thin dashed 
lines forming the light blue and green error ranges. To indicate 
the resulting range of maximum grain sizes, we label the cross¬ 
ing points with the observed values and their error range derived 
from the off-measurement uncertainty by circles. Bars near the 
x-axis aid in reading the ( 22 -range corresponding to the agree¬ 
ment of observed and theoretical coreshine. As is visible in the 
figure, the bars partially overlap, indicating that by extending the 
distribution up to sizes around 0.84 jum, the model reproduces 
the observed data. 

I Andersen et al.1 (l2013l) have modeled the 2.2, 3.6, and 8 jum 
SFB along a cut through L260, and found maximum grain sizes 
around 1 jum, which is close to what we find for our ice-coated 
silicate and carbonaceous grains. Their calculations required 
considering an enhanced incident DIRBE field boosted by a fac¬ 
tor of 1.7, which might have been due to an error in the older RT 
code and is no longer needed in the new code version. 

In Panel B, we show the ratio of the blue and green curves 
from panel A, which is the theoretical band ratio R with the 
same notations, and the bullets indicate the crossing points of the 
model with the observed maximum coreshine SFBs. We show 
just the bullets and not the dashed lines from Panel A to be able 
to plot more than one R curve in the following figures without 
confusing the plots. As expected, the curves show the increase 
due to the pole in R when the 4.5 jum SFB cannot overcome the 


^ Diffuse Infrared Background Experiment, see 
www.lambda.gsfc.nasa.gov/product/cobe/dirbe_overview.cfm 
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maximal grain size a2 [iiim] 



maximal grain size a2 [iiim] 

Fig. 4. Panel A: Background-subtracted theoretical surface 
brightness through the core center at 3.6 (4.5)//m in blue (green), 
respectively (see Eq. 13 as a function of the maximum grain size 
of the dust size distribution for the core L260. The roots of both 
functions indicate the regions where the core becomes visible 
in excess of emission for each wavelength (region outside yel¬ 
low regions). The dashed lines give the observed maximum core¬ 
shine value with its error range, and the dots label the crossing 
points with the theoretical curve. Panel B: Resulting theoreti¬ 
cal ratio Rtheo of the surface brightnesses shown in Panel A as 
a function of the maximum grain size. The dots from Panel A 
mark the agreement of the central surface brightnesses with the 
observed range for each wavelength. 


background extinction, and the core turns to an extinction ap¬ 
pearance at smaller maximum sizes. 

In the optically thin case, the curve in Panel B should be in¬ 
dependent of assumptions about the spatial density distribution 
within the core, while the SFB profiles in Panel A assume a cen¬ 
tral column density or optical depth and therefore rely on spatial 
parameters. For optical depths higher than unity both panels will 
be affected by the assumptions made within the density distribu¬ 
tion model. 

We note that R also depends on the position of the core since 
the ISRF is anisotropic with a maximum at the GC, and con¬ 
volving it with the anisotropic phase function will lead to dif¬ 


ferent values for different core positions. As a technical note we 
add that when referring to the theoretical quantities like the opti¬ 
cal depth, we skip the notation, indicating the integral over 
grain size to improve readability from this point on. 

4. Coreshine modeling of the core sample 

The general approach taken in this section is to compare the mea¬ 
sured R distribution in a region of the core with the theoretical 
ratio Rtheo derived from a model core. The core properties and 
the maximum grain size are varied without exceeding the limit¬ 
ing optical depth (Sect. iTTl) . Since two sets of 3.6 and 4.5 yum 
images can have the same ratio map but might differ strongly at 
each wavelength, we also verify the agreement of the SFBs. For 
cores with r < 2, it is sufficient to consider a single SFB value. 
We have chosen the maximum observed SFB and not some aver¬ 
age for two reasons. First averaging considers more regions with 
worse S/N, increasing the overall error, and second the averaging 
would bring in more assumptions about the way the averaging is 
done. Even if the core has regions with r > 2, we try to select 
a region in the observed and model image that has optical depth 
below that limit. In this way we take advantage of the optically 
thin approximation being valid, which removes spatial depen¬ 
dences from the problem. For cores with r > 2, we perform ’’full 
RT” modeling by calculating images, building ratio maps and 
an R distribution, and fitting the observed R distribution and the 
central column density within the optical depth limit. 

Throughout the section, we make use of the approximate 
core properties summarized in Table [T] 

4.1. Modeling details explained: the core L260 

To describe the modeling approach in detail, we have chosen the 
co re with the brightest core shine in the sample: F260. According 
to ISteinacker et al.l (l2014ah . the core is located in a region fa¬ 
vorable to coreshine detection because it is near the GC in the 
PoSky, but the latitude is high enough to prevent a strong back¬ 
ground. Still it has the second strongest 3.6 jum background of 
the cores modeled in this paper. 

Coreshine is a weak extended feature, and especially the 4.5 
jLim SFB is only factors 2-8 larger than the sum of instrumen¬ 
tal and point source noises. Correspondingly, we have selected 
special regions in the images to minimize the impact of noise 
and stellar point spread function (PSF) on the analysis. The two 
background SFBs and entering Rtheo are taken from the 

DIRBE maps by subtracting the stellar flux measured by WIS^. 
We did not subtract a diffuse foreground component of the ra¬ 
diation field used as background of the core sin ce the current 
mode ls give only indications (for a discussion see lFefevre et al.l 
l2014l) . We briefly address its effect in Sect. [5] For 3.6 pm, the 
bands of IRAC and DIRBE are comp arable, for 4.5 jum we inter¬ 
polate from the other bands (see also I Andersen et al.]l2Q13l) . The 
sum of the background and foreground contributions can be de¬ 
termined from the SFB near but off the core (”ofF’ characterizes 
quantities derived in this way) 

z) = Itp, z) + z). (13) 

The details of how this subtraction is carried out are explained in 
Appendix We note that the subtraction removes the possible 
zodiacal light contribution. 


3 www.nasa.gov/mission_pagesAVISE/ 
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Fig. 5. Data analysis for L260. Warm Spitzer IRAC surface brightness maps (off-subtracted) at 3.6 and 4.5 pm (Panels A,B) and 
cold Spitzer map at 8 pm (Panel C) in a 5.5'x5.5' region around the core. The white square indicates the l.l'xl.l'-region where 
surface brightness ratios have been measured. The optical depth for extinction of external radiation passing through the core and to 
the observer (directionally averaged) for a core with the overall properties of L260 is shown in Panels D and E (3.6 and 4.5 pm). 
Panel F: Horizontal (thick) and vertical (thin) cuts through the central 200 pixel in the two bands (black for 3.6 pm and green for 
4.5 pm). Panel G: Map of the ratio of the off-subtracted surface brightnesses R at 3.6 and 4.5 pm. Panel H shows the R pixel number 
distribution within the white square (pixels affected by point sources or by the column-pulldown effect have been masked). Panel I: 
Rtheo as a function of the maximum grain size a 2 (thick black), with the notation from Fig.|4l The observational R-distribution from 
Panel H is plotted and color-coded, and both the maximum value and the two half-maximum values are plotted as dashed lines. The 
horizontal black bar indicates the a 2 agreement range of Rtheo and Robs^ 


Figure [5] shows the steps in the analysis for the core F260. 
The off-subtracted SFB maps at 3.6 and 4.5 pm in a 5.5'x5.5' 
region around the core are shown in the Panels A and B, respec¬ 
tively. The cold Spitzer 8 pm map is given in Panel C to indi¬ 
cate the location and the shape of the core (”Coreplanets_cores”: 
Program ID 139, PI N. Evans). L260 is seen in emission in 3.6 
and 4.5 pm. The white square indicates the region where SFB 
ratios are measured. Since the core shows no sign of central de¬ 
pression in the two bands, the square is placed in the central 
region of extinction in the 8 pm map. The size and precise po¬ 
sition of the square have been adopted to avoid the impact from 
nearby stars and their PSF, to minimize the impact from column 
pull-down in the image, and to recover the grain information in 
the part of the core where the coreshine signal-to-noise is max¬ 
imum in both bands. Alternative shapes of the measuring area 


would have been possible, but an investigation of the ratio maps 
(Panel G, see farther down) showed that the results are not sen¬ 
sitive to the shape of the area. All remaining pixels that show 
impact from a stellar PSF or column pull-down are masked and 
not modeled. 

To check the distribution of optical depths for extinction in 
the region where we measure R, we show the direction-averaged 
optical depth map of a core with the overall properties of F260 
in Panels D and E. The direction average takes into account that 
radiation is extincted while reaching the FoS from outside and 
on its way to the observer. As simple core model, we use an el¬ 
lipsoidal core with two of the three axes identical and a radial 
density profile that flattens in the inner part while following a 
power law with index -1.5 outside. Core mass and outer radius 
are taken from TablelU and the kink radius is assumed to be one- 
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third of the outer radius. We note that both the density power-law 
index and the kink radius will depend on the evolutionary stage 
of the core, but more evolved cores might have a smaller kink 
radius, for example, and for LI544 we have used a different ra¬ 
dial profile. For L260, the mean optical depth values reach about 
0.4 for 3.6 yum. We there fore use single-scattering RT to model 
the core. We note that the models shown in Panels D and E for 
the core sample do not aim to model the spatial density distribu¬ 
tions but only serve to estimate the optical depth variation, with 
the exception of the cores with maximum optical depth above 2 
where we create images in full RT based on the shown spatial 
model. 

To read approximate values of the SFB band ratio, the noise 
level and the contamination by stars, we show in Panel F hori¬ 
zontal (thick) and vertical (thin) cuts through the center of the 
white line frame for 3.6 (4.5) jim with black (green) lines. While 
there are local gradients along the cuts that vary the surface 
brightness ratio, the figure already allows reading an approxi¬ 
mate ratio of about 2 to 4. The cuts also indicate the level of 
stellar radiation impact and the overall noise of of about 0.005 
MJy/sr. Figure |4| shows the two maximum surface brightness 
values and their error bars as dashed vertical lines. The error 
bar contains not only the error read from Panels A, B, and F in 
Fig.B but also the error in the off-measured surface brightness 
that is subtracted from the measured surface brightness. The en¬ 
tire map of off-subtracted surface brightness ratios R is shown 
in panel G. Most variations in the white square are due to stellar 
sources. We masked all pixels affected by stellar PSFs and by 
column pull-down. In panel H, we show the R pixel number dis¬ 
tribution across all unmasked pixels in the white square. It has a 
maximum at R =2.37, and as a measure of its size, we use the 
approximate width of the distribution at the half maximum value 
indicated by thin vertical lines (the full width at half maximum, 
FWHM, is about 1). 

In Panel I, the range of observational ratios Robs is compared 
to the theoretical ratio Rtheo as a function of a 2 as derived in 
Sect. [3l with the background values given in Table [2] and us¬ 
ing the core sky position for the phase function integration of 
the scattering integral. The limited optical depth approach is ap¬ 
plied, which prevents the maximum optical depth for extinction 
to exceed 1 at T = 2.2 //m (see Table[T]). The green color-coding 
indicates the observed R-distribution from Panel H. 

The crossing points of the theoretical R curve with the ob¬ 
served coreshine range at 3.6 yum (blue) and 4.5 yum (green) are 
shown as dots as in Fig. |4] with small dots indicating the error 
range and large dots giving the crossing point with the exact 
value. There is an overlap of the ^2 ranges of the two bands in¬ 
dicated as blue and green horizontal bars at the bottom of the 
figure and the black bar marking the < 22 -range where Robs and 
Rtheo agree. Thus the model is able to meet all three observational 
constraints for maximum grain sizes between 0.78 and 1.02 yum. 

We have performed extensive tests on varying the 
spectral slope of the siz e distribution as suggested by 
IWeingartner & Drain3 (1200ll) . The lower grain size limit has no 
effect on the coreshine si nce the scattering is dom inated by the 
largest grains as shown in ISteinacker et al.l (l2010l) . The effect of 
the size distribution slope on the derived grain size limits is as 
expected: flattening the distribution will decrease the limits, as 
the number of bigger grains increases with a flatter slope. Since 
they are responsible for most of the observed coreshine, less big 
grains are needed to meet the observed flux. The impact on the 
limits is moderate since all opacities are integrated over the size 
distribution. For a slope of -2 instead of -3.5, the limits move by 
about 0.15 yum. More importantly, the overlap between the ^ 2 - 


ranges derived from the two observational constraints remains 
the same. As a result, we argue that the derived grain size limits 
may vary by several tens of percent when changing the slope of 
the size distribution power law in a reasonable range of [-2 to -5], 
but this will not affect the conclusions. For the models presented 
in this section, we have used a spectral slope of -3.5. 

4.2. ecc806 (PLCKECC G303.09-16.04) 

The second brightest core in the sample, PLCKECC G303.09- 
16.04 (we will use the shorter name ecc806 here), is surrounded 
by more filaments with coreshine that are located outside the 
area shown in Fig. [6l As for L260, the position favors the de¬ 
tection of coreshine because it is close to the GC in longitude 
and far below (b=-16.04°) the Galactic plane. Moreover, it has 
the lowest background in the sample in both bands. There is no 
sign of depression in the images and in the cut in Panel F. We 
therefore use single-scattering RT. The chosen white frame en¬ 
closes the core center with optical depths below the r limit. In 
Panel C, the 250 yum map observed by Herschel (observation 
id 1342213209) as part of the Gould Belt Survey (PI: Andre) 
is shown for comparison. As with the extinction pattern at 8 
yum shown for the other sources, the thermal emission pattern 
also matches nicely the coreshine pattern. The SFB profiles are 
remarkably flat in the central part. The R distribution is only 
slightly broader than that of L260 (FWHM about 1.2) but with a 
maximum at larger R near 3. The fit indicates that the maximum 
grain size in ecc806 is around from 0.4 to 0.52 yum. 

4.3. L1262 

Modeling of the starless core in the binary core L1262 (CB 244) 
is complex owing to the nearby YSO with its outflows creating 
H2 emission in the shocked gas in both bands. L1262 also has a 
remarkably extended cloudshine that is visible at shorter wave¬ 
length and is the most distant core in the sample (about 200 pc). 
The YSO and its warm-dust surroundings may give rise to an 
additional radiation field. Since the density distribution is com¬ 
plex, it is hard to model this component in order to derive the 
intensity and directional dependency at the location of the star¬ 
less core. With 60° away from the GC, the phase function of the 
largest grains will not increase coreshine by beaming the strong 
GC radiationforward. 

The core position in the PoSky is elevated enough (Z?=12.4°) 
to have a low background, especially at 4.5 yum, which allows the 
detection of the weak coreshine signal in this band. The mass of 
the core is rather uncertain with a mean value of 5 M© and an 
error bar of 2 M©. However, the column density range given in 
Table [T] also points to r > 2, and the 3.6 yum image in Panel 
A reveals a depression right where the extinction against the 
background has a maximum in the cold 8 yum image in Panel 
C (’’darkcloudsJRAC”: Program ID 94, PI C. Lawrence). We 
therefore have chosen maximum optical depth values around 4.5 
(scaled with the mean column density values of the other cores) 
and use full RT. We have chosen a white frame away from the 
YSO and its extended tail and excluded the part of the square 
where the extinction seen in Panel C of Fig. [7] indicates the cen¬ 
ter of the core. As mentioned before, this improves the modeling 
because at lower optical depths, some of the spatial dependen¬ 
cies cancel out. The observed R distribution is about a factor of 2 
narrower than for L260 with a lower mean value of 1.75. Aside 
from the color-coded observational R-distribution, Panel I shows 
the distribution of the theoretical R-values with a thick line for 
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Fig. 6. Data analysis for ecc806 (PLCKECC G303.09-16.04). For a detailed explanations of the panels see Fig.[5l Panel C shows 
the SFB at 250 jjm as observed by Herschel. 
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the mean value and thin lines for the FWHM. The a 2 ranges 
overlap between 0.52 and 0.71 yum. 

4.4. L1517A 

The FI517 cl oud in Taurus contains four filaments with five em¬ 
bedded cores (iHacar & Tafallall2QTl1) . For the modeling we have 
chosen core A since it was visible in both bands in our new deep 
observation images. The coreshine of all Taurus cores (in our 
sample F1517A and F1512, F1544, F1506C, and F1498 with 
I around 180° discussed later) benefit from a side maximum 
in the dust- model phase function of large grains that predicts 
enhanced backward sc attering of the strong GC radiation (see 
ISteinacker et al.|[^14ah . 

The core is located at the border of the 3.6 jim image (Fig.O 
see also Panel C, ’’darkcloudsJRAC”: Program ID 94, PI C. 
Fawrence), but were able to select a region that includes the core 
center and most of the core. We considered using cold 3.6 jum 
data instead, but decided to stay with the warm data for the com¬ 
parison to minimize the impact of systematic effects between 
the observing runs because the data show no strong increase in 
noise compared to the warm 4.5 yum data. There is no indication 


of any depression, and the mass and column densities point to¬ 
ward a maximum optical depth value for extinction below 1. We 
therefore use single-scattering RT. The cuts show a clear con¬ 
stant factor between the two bands in both cut directions, as one 
would expect from optically thin scattered light where the SFB 
is proportional to the column density. This should give rise to a 
narrow R distribution. On the contrary. Panel H shows a broader 
R distribution than F260 (FWHM is about 1.5, peak at 2.75). 
We attribute this to the source having the strongest background 
in the sample. The distribution also is asymmetric with a strong 
tail toward higher R. Panel G shows an increase in R values in 
the upper part of the white frame Probably because of the loss 
of the 4.5 yum coreshine signal, and we attribute the asymme¬ 
try in Panel H to this signal loss. The model is able to meet all 
observational constraints for a 2 values above 1.1 yum. 


4.5. L1512 

This Taurus core is located close to the Galactic plane and has 
an enhanced background (second strongest 3.6 yum background 
in the sample). As for FI517A, we use the new warm data show¬ 
ing the core at the image border instead of using center-image 
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Fig. 7. Data analysis for LI262. For s detailed legend see Fig.[5l The binary core harbors a YSO to the left. The core was located 
at the IRAC 1 image border leading to a dark zone in Panels A and G. (The R values in the dark right zone are not meaningful.) 
Since the modeling is performed using full RT, Panel I compares the observed R distribution as color-coded and the theoretical 
R-distribution with a thick solid line for the mean value of the theoretical R distribution and half maximum values as thin lines. 


3.6 fim cold Spitzer data in Fig.O The cold Spitzer map at 8 pm 
shown in Panel C is from the ’’darkclouds JRAC” program (ID 
94, PI C. Lawrence). Based on the core property estimates from 
Table [T] the optical depth analysis indicates that the central part 
of an ellipsoidal representative of LI512 would see optical depth 
above 1, and therefore extinction and second scattering make full 
RT necessary. The image in Panel A and cuts like the one given 
in Panel F show little central depression though, but strong stel¬ 
lar PSFs make this analysis difficult. We therefore performed full 
RT modeling. The 4.5 pm signal is just about a factor 2 above the 
noise with the image showing just hints of excess of emission. 
The R distributions are correspondingly broad (FWHM of about 

2.6 around R=3.6) and conclusions from the modeling must be 
made with caution. The full RT distribution for R overlaps with 
the observed R values for maximum grains between 0.48 and 
0.61 pm. 

4.6. L1544 

The Taurus core LI544 is one of the few known cores across 
the sky that show a prominent inner depression of SFB in both 


bands. While the mass seems to be comparable to that of LI512 
showing no clear depression, the central SFB of LI544 very 
likely drops because of optical depth effects of the centrally con¬ 
densed radial profile. The 8 pm map shows a strong extinction 
pattern at the core center (”cores2deeper”: Program ID 20386, 
PI P. Myers). This would mean that we can expect moderate op¬ 
tical depth in the outer parts of the elliptically shaped core and 
a strong rise in r near the center. Using a profile with a power- 
law index of -2.5 in the outer parts and a flattening starting at 
the kink radius of 1 kau, the mean optical depth maps shown in 
Panels D and E in Fig. [TOl pro vide such a pattern. Scaling r with 
the estimated central column density of LI544, we use a r limit 
of 6 at 4=3.6 pm. The white square region was placed outside 
the depression, and we use full RT. The derived ratios are low 
with a mean value of 1.25, and the FWHM is comparable to that 
of L260 (about 0.95). But as is visible in Panel F, the data are 
extremely noisy with a dense stellar background field interfer¬ 
ing and with a loss of the 4.5 pm signal in the lower white frame 
part that causes a wing in the R distribution at high R. The model 
suggests large grains beyond a2il.5 pm to reproduce the SFBs 
and ratios. 
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Fig. 8. Data analysis for LI517A. For detailed legend see Fig.[5l The core was located at the IRAC 1 and 2 image borders leading 
to a dark zones in Panels A,B, and G, (Here the R values in the dark lower and right zones are not meaningful.) 





4.7. L1506C 

This extended Taurus core with low dens ity, low turbulence, 
but strong depletion (iPagani et al.l l2010ah was discussed in 
ISteinacker et all (l2014bh as a core that either is s till contract¬ 
ing to become a prestellar core (as suggested by iPagani et~n] 
l2010ah or has passed through a core phase and is dissolving 
again. Since the core is invisible at 8 yum, w e show in Fig M in¬ 
stead a MAMBO II map (reproduced from IPagani et al.|[2010ah 
at 1.2 mm revealing the same pattern (in thermal emission) as the 
coreshine. The large Galactic latitude of -17.57° helps to make 
the coreshine visible. The 4.5 yum signal is just a factor 2 over the 
noise but the entire core pattern is nicely visible and exhibits no 
strong impact of potential noise confusion. From the noisy cuts 
displayed in Panel F, it is difficult to judge on potential central 
depressions, since the images in Panels A and B do not reveal a 
darkening, and the basic parameters suggest optical depths be¬ 
low 1 in the center. As a result, we placed the white frame near 
the core center, thereby avoiding the two large stellar PSFs, and 
use single-scattering RT. The R map in Panel G shows the inter¬ 
esting feature that the entire clump visible in the image seems 
to have R values comparable to the mean value measured in the 
white frame as 1.7 (FWHM about 1). 


The model is able to account for the SFBes based on grains 
with maximum sizes around 0.4 yum, while the ratios produced 
by such grains of about 7 would be too high to match the 
observed low values between R=l.2 and 2.2. This steep ratio 
points toward grains wit h sizes larger than 0.7 yum . The analysis 
of L1506C presented in ISteinacker et al.l (l2014bl) using a grain 
growth model and simple spatial mixing of the size distribution 
pointed toward large grains beyond 1 yum in size. However, in¬ 
cluding the spatial mixing in the growth process would likely 
reduce the efficiency of the growth process, leading to smaller 
maximal sizes. 


4.8. L1439 

The core L1439 (CB26) is located only 6^ above the Galactic 
plane in the PoSky and is therefore seen against a crowded star 
field. The 4.5 yum signal is only a factor 2 above noise (Fig.fT^. 
At its rim, a YSO creates a jet, and it is not clear whether the 
stellar component or the warm dust adds local radiation to the 
ISRF. The cold Spitzer map at 8 pm is shown in Panel C (’’dark- 
clouds JRAC”: Program ID 94, PI C. Lawrence). The core pa¬ 
rameters listed in Table \T\ make central optical depths below 1 
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Fig. 9. Data analysis for L1512 (CB27). For detailed legend see Fig. [5] Panel A contains an image border that leads to mean¬ 
ingless ratios in Panel G in the dark righthand area. Panel I compares the observed (dashed) and theoretical (solid) R distribution 
characterized by the mean value (thick) and the FWHM values (thin). 
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likely. We therefore placed the white frame right on the core and 
used single-scattering RT. The modeling with a variable r limit 
led to values around 0.5 for 3.6 jum. The core is located at the 
border of the 3.6 yum frame so that the dark righthand part of the 
R map displayed in Panel G is meaningless. The R pixel number 
distribution shown in Panel H of Fig. [T2| contains an asymmet¬ 
ric wing from the map parts with a weak 4.5 yum signal. It has a 
maximum around 1.7 and a FWHM of 1.4. Panel I reveals that 
the model produces no common size range to meet the observed 
SFBs and the observed low Robs ratios. 


4.9. L1498 

LI49 8 in Taurus is di scussed to be on the verge of becoming a 
core (iLada et al.ll2004l) . and only a faint extinction pattern is vis¬ 
ible at 8 yum in Panel C of Fig.[T3](”darkclouds JRAC”: Program 
ID 94, PI C. Lawrence). Positioning of the white frame was 
guided by the 3.6 yum image border so that it was placed in one 
half of the ellipsoidal shape of the core. The R maps in Panel 
G show that this choice does not largely affect the results since 
the other parts of the core seem to have similar values to those 


of the core L1506C. The core parameters suggest central opti¬ 
cal depths in the range of L1506C, so we use the same optical 
depth limit and single-scattering RT. Given the low 3.6 yum sig¬ 
nal, LI498 shows a surprisingly strong 4.5 yum SFB yielding to 
low ratios around 2 (FWHM is 1.55). The model grains produce 
too much coreshine when using the grains with size above 0.77 
yum as suggested by the model to fit the Robs distribution so that 
the a 2 ranges do not overlap. 

4.10. RhoOph 9 

The core Rho Oph 9 is located in the most complex environment 
of all cores in the sam ple (Gal, coord. 354.3 7 -f 16.17). The cold 
Spitzer data shown in iPagani et"^] (1201 Obh show nearby PAH 
regions in all four bands with reduced SFB at 4.5 pm owing to 
the lack of a main PAH feature in this band. 

The excess of emission that is seen near the strong extinction 
pattern is also present at 4.5 pm similar to the pattern seen in 
LI544. Unlike the PAH emission, it disappears for larger wave¬ 
length in the cold Spitzer data as expected from scattered light. 
Figure [14] shows the 3.6 and 4.5 pm warm Spitzer images. The 
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Fig. 10. Data analysis for L1544. For detailed legend see Fig.[5l The core shows an inner surface brightness depression, and the 
white frame was chosen outside this region. The surface brightness signal is noisy owing to the many stars in the field. 
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overall pattern is similar to LI544 with a strong central SFB 
depression and an ellipsoidal emission region around. But dis¬ 
entangling the coreshine emission from variations in the back¬ 
ground and PAH emission turned out to be difficult. We have not 
modeled the core since we could not determine an off-region that 
was clearly free of PAH emission and near the core, and hope to 
return to the source when more data are available. We show the 
maps and the ratio map for completeness. 

The peak and background SFBes of the modeled sources are 
summarized in Table [2] along with the derived observed SFB ra¬ 
tios Robs^ 

5. Discussion 

In Fig. [15] we summarize the a 2 ranges derived from the model¬ 
ing with the same notation as in Fig. [5j blue, green, and black 
bars for agreement with the 3.6 yum maximum SFB, with the 4.5 
Jim maximum SFB, and with the observed SFB ratio Robs, re¬ 
spectively. 

For L260, ecc806, L1262, L1517A, L1512, and L1544, the 
bars overlap so that the model is able to reproduce the observa¬ 
tions with a single grain size distribution. The derived maximum 
grain sizes for these sources are around 0.9, 0.5, 0.65, 1.5, 0.6, 


and 1.5 yum, respectively. The value derived fo r L260 of about 
0.9 yum is consistent with the results derived in I Andersen et al.l 
(I2ni3h . 

The grain size limit of about 0.5 yum in ecc806 is surprisingly 
low given that it is bright in coreshine (about 0.075 MJy/sr) and 
that large grains dominate the scattering. At a longitude of about 
303°, radiation from the GC already needs large scattering an¬ 
gles to reach us. Because grains larger than 0.5 yum scatter mainly 
in the forward direction, their contribution to the observed signal 
is reduced. However, the low background due to its high eleva¬ 
tion (b about -16^) aids for seeing a strong coreshine. 

The structure around the starless core in the binary core sys¬ 
tem LI262 is complex and not represented well by an ellipsoid. 
Correspondingly, the derived images carry the errors from the 
spatial density because of both shading and multiple scattering. 
Nevertheless, maximum grain sizes of 0.65 yum agree with the 
measured ratios and maximum SFBs in the chosen region out¬ 
side the core center. They may actually be smaller since the 
nearby YSO could contribute to the local radiation field both by 
direct irradiation and by causing hot dust and PAH emission in 
its vicinity. Further exploration with a better spatial model may 
answer that question. 
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Fig. 11. Data analysis for L1506C. Panel C shows the MAMBO II map reproduced from iPagani et alJ (l20I0ah . For detailed legend 
see Fig. [21 


Table 2. Peak and background surface brightnesses at 3.6 and 4.5 jim and derived observed coreshine ratios R. 


Source 

j3.6 

peak 

[MJy/sr] 

[MJy/sr] 

t3.6 

hg 

[MJy/sr] 

C 

[MJy/sr] 

R 0 Z 75 

L260 

0.085 

0.04 

0.098 

0.114 

1.9-2.9 

Ecc806 

0.075 

0.022 

0.026 

0.036 

2.3-3.6 

L1262 (CB244) 

0.070 

0.035 

0.031 

0.042 

1.4-2.2 

L1517A 

0.055 

0.02 

0.150 

0.171 

2.0-3.7 

L1512(CB28) 

0.035 

0.01 

0.052 

0.082 

2.5-5.0 

L1544 

0.030 

0.02 

0.040 

0.090 

0 . 8 - 1.8 

L1506C 

0.027 

0.015 

0.032 

0.086 

1 . 2 - 2.2 

L1439 (CB26) 

0.025 

0.01 

0.041 

0.066 

1.1-2.4 

L1498 

0.025 

0.01 

0.046 

0.090 

1.3-2.9 


Strong coreshine (about 0.055 MJy/sr) also comes from 
LI517A but on top of a much stronger background. 
Correspondingly the grain size limit a 2 i 1.5 yum derived by the 
model is larger than, for example, that of ecc806 with similar 
maximum SFB values to enable a stronger scattering. Like all 
cores in Taurus, L1517A benefits from side maxima of the phase 
function supporting back scattering of the strong GC signal. 

For the cores with the lowest coreshine SFB L1506C, L1439, 
and LI498, the three different a 2 ranges do not overlap in one co¬ 


herent way. The maximum grain sizes to explain the ratios are in 
all cases larger with just a lower limit. To investigate whether 
this systematic trend arises from a simplification of the model, 
we consider a local radiation field in addition to the ISRF. Since 
the cores are located in star-forming regions where warm dust 
is present near the YSOs, such a local radiation field is not un¬ 
likely. Moreover, the DIRBE map reveals radiation on top of the 
stellar component from the regions where cores are located, with 
a spectrum increasing with wavelength in the considered bands. 
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Fig. 12. Data analysis for L1439 (CB26). The 3.6 jim image (panel A) contains an edge which creates a dark region in the ratio 
map in Panel G where the ratio is not defined. For detailed legend see Fig. [S] 
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While the radiation sources near the core will also contribute 
to the DIRBE fiux from that region, if it is not shielded along 
the LoS, sources very close to the core will make an increased 
contribution to the local field due to the squared distance depen¬ 
dency of the fiux. 

The directional dependence of the local field is unknown and 
is likely to be different from that of the ISRF. To explore the im¬ 
pact of a local field, we assume a simple field with no directional 
variation and an intensity of 0.05 MJy/sr (about 1/6 of the mean 
DIRBE field) at 3.6 yum. We modify the spectral shape by a fac¬ 
tor (4/2.2 [imY so that a variation in 6 takes the rise in the local 
field intensity with wavelength into account. Figure [T6| shows 
three cases for the core LI439. The thin line gives the theoreti¬ 
cal R without a local field. When adding a constant field without 
spectral variation {6 = 0, medium thick line), more radiation is 
scattered toward the observer. The grains can thus be smaller to 
explain the observed SFB, but the overlap between the ranges 
for each band, and the ratio is not improved. Assuming a local 
field with an increase described by ^ = 1.5, this trend continues to 
lead to no improvement in the fit. For cores with an YSO in the 
vicinity or even embedded, a local anisotropic radiation field is 
added to the incid ent field, and i t s impa ct on the coreshine ratios 
was discussed in iLefevre et al.l (120141) . They find higher core¬ 


shine SFBes and lower observed coreshine SFB ratios h.elh.s 
for embedded cores than for starless cores. This would make it 
more difficult in our model to converge to a common maximum 
grain size. A lternatively, the s ize dist ributions explored in that 
paper and in ISteinacker et al.l (l2014bh mark the beginning of a 
more systematic search of the large size distribution parameter 
space for a dust model that can explain the presence or absence 
of coreshine i n the various ob s ervatio nal bands. 

For L260, lAndersen et al.l (l2013l) also used the observed Ks 
band to constrain the size distribution. At shorter wavelengths, 
the central depression is more likely. It needs to be verified that 
peaks in the grain size distribution at sizes around 1 yum will 
be compatible with the coreshine and extinction observed at 
shorter wavelengths. This emphasizes the importance of deep Ks 
band and NIR observations in general of cores with coreshine. 
Moreover. lLefevre et al.l (l20l4 argue that both the NIR/MIR ra¬ 
tio for the core outer layer and the absence of emission at 5.8 yum 
for any layers eliminate a mix of silicates and carbonates that 
both include grains above about 2 yum in meaningful quantity. 

A word of caution is needed for the uncertainties in the back¬ 
ground estimates. A better estimate of the foreground diffuse 
radiation field would improve the current uncertainty that en¬ 
ters coreshine modeling. Weaker background would increase the 
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Fig. 13. Data analysis for L1498. For detailed legend see Fig.[5l 



Fig. 14. 3.6 and 4.5 yum maps, and ratio map for Rho Oph 9. 


coreshine SFB in both bands so that smaller grains would be suf¬ 
ficient. The cores where the model fails would req uire the oppo¬ 
site. M oreover, the background values derived by iLefevre et al.l 
(l2Q14l) using a different method can vary from one region to the 
other by a factor of up to two. This can modify the results for the 
cores with weak coreshine or strong background. For future 3D 
modeling of a single source, it could be important to put more ef¬ 
fort into modeling the foreground and background contributions. 


One particular finding about L1506C and LI498 is interest¬ 
ing. Both cores show little variation in Robs across the entire core 
as far as 4.5 yum coreshine is detectable. For both cores this ra¬ 
tio is small, as expected from the presence of larger grains. Both 
are discussed as being on the verge of becoming low-mass cores 
with typical properties, and also show central gas depletion that 
requires a density that might not be reached in these cores. 
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Fig. 15. Summary of the derived maximum grain size ranges for 
the modeled sources. The black bar gives the range where Rtheo 
is in the FWHM range of Robs^ and the blue (green) bar shows 
the range for the model maximum coreshine to be within the ob¬ 
served maximum coreshine range at 3.6 (4.5) yum, respectively. 
For L1439, we also give size ranges for a model that includes a 
local radiation field. Light green boxes indicate ranges of agree¬ 
ment. 
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Fig. 16. Impact of a local radiation field with no directional vari¬ 
ation and a wavelength dependency following {A/l.ljdmY on the 
coreshine ratio of LI439 as a function of the maximum grain 
size. 


For both cores the existence of larger grains would be sur¬ 
prising if the grains are to be formed by coagulation since their 
densities are not high enough to co agulate them, as indicated by 
studies of ISteinacker et al.l (i2014bh for L1506C. One of the hy¬ 
potheses presented in that work was that L1506C went through 
a phase of higher density and stronger turbulence to explain the 
existence of larger grains and strong gas depletion. With the sim¬ 
ilarities between the two cores also emerging from the work pre¬ 
sented here, this could also be the case for LI498. 


6. Summary 

This work is based on deep 3.6 and 4.5 yum IRAC Spitzer warm 
mission data taken during Cycle 9. The sample is composed of 
ten cores, nine of which show coreshine at both 3.6 and 4.5 pm, 
and one where coreshine is unconfirmed owing to strong PAH 
emission. The observations and the work in this paper were mo¬ 
tivated by the expectation that coreshine observed at MIR wave¬ 
lengths might be dominated by dust grains in the micron size 
range and that opacity models predict a stronger impact of the 
largest grains at longer wavelengths. 

Given the vast parameter space of the general grain scatter¬ 
ing problem in cores caused mainly by the unknown density 
structure, the presented model approach makes use of the fact 
that in the optically thin limit, both the observed coreshine in¬ 
tensity and the intensity losses by extinction through the core 
are proportional to the column density of grains. By considering 
ratio maps of 3.6 and 4.5 pm intensities, the impact of the spatial 
density distribution is expected to be reduced for the cores with 
moderate optical depth < 1. Moreover, we use the constraints 
from the observed SFBes based on the expected column densi¬ 
ties. Since the core properties are known only within factors of 
a few, and the opacity changes caused by increasing the grain 
size can be substantial, it is essential for the modeling not to 
exceed the observed optical depth limits, which are provided in 
the form of depression of the central SFB for cores with opti¬ 
cal depth for extinction above about 2. When the column density 
is varied within its uncertainties, we make sure that this optical 
depth limit is never exceeded. 

We applied a model based on the opacities (cross sections 
and phase functions) of ice-coated silicate and carbonaceous 
spheres and a size-distribution following a power law (index - 
3.5) ranging from 0.01 to a variable maximum grain size for both 
species. The incident field is approximated by DIRBE maps, and 
the radiation right behind the core is constructed from the ab¬ 
solutely calibrated DIRBE map with point source contribution 
subtracted using the WISE catalogue. 

We outlined three ways to perform RT and describe the gain 
and potential errors by comparing them. We illustrated that the 
spatial and opacity information is blended into the derived SEB 
by the action of multiple scattering and extinction for increasing 
optical depth. Eor an example core, we found that the error of 
the optically thin approximation is about 14% compared to full 
RT as long as the optical depth for extinction stays below 1. The 
theoretical coreshine ratios are described and applied for a model 
core with the properties of L260. They have a singularity for the 
limiting grain size where the net effect of core extinction and 
scattering is zero and then decline toward -values around 1 
for grains beyond 2 pm in size. 

We described the modeling procedure in detail for L260, 
also using longer wavelength maps like cold Spitzer data at 8 
pm to construct a simple spatial model. The model was used 
to estimate the optical depth variation. Selecting a region in the 
ratio map of the two bands that primarily shows optical depth 
smaller 1, the pixel distribution as a function of the ratio serves to 
characterize the observed ratios. Then we applied either single¬ 
scattering or full RT depending on the assumed maximum opti¬ 
cal depth to derive theoretical R values or distributions to com¬ 
pare them to Rtheo along with the constraints from the observed 
maximum coreshine SEBes. 

Eor the cores L260, ecc806, L1262, L1517A, L1512, and 
LI544, the model is able to account for the observed central 
SEBes and the ratio maps with maximum grains size around 0.9, 
0.5, 0.65,1.5, 0.6, and ^1.5 pm, respectively. The low grain size 
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limit for ecc806 is remarkable given the bright coreshine from 
this core, and might be related to the low background. 

The maximum grain sizes found agree with the findings of 
earlier studies listed in the introduction, suggesting that core¬ 
shine indicates the presence of grains larger than found in the 
diffuse ISM. No obvious correlation with basic core properties 
is evident for our sample of nine modeled cores. 

The model fails to simultaneously reproduce the SFBes and 
ratio maps of L1506C, L1439, and L1498. The grain size limits 
derived from the ratio maps are larger than the limits from the 
coreshine intensity. Since the coupling to the spatial structure is 
moderate (Text < 1), this possibly points to assumptions that are 
too simplified or to the action of an additional component. 

Assuming a constant local field with a rising spectral shape 
on top of the ISRF described by the DIRBE all-sky map does 
not improve the modeling for LI439. The fiat R distribution of 
the cores L1506C and LI498 gives rise to speculations that both 
cores host larger grains across the core, either as a primordial 
component or created in a former more dense a nd more turbulent 
phase as discussed in ISteinacker et al.l (l2014bh for L1506C. 

Rho Oph 9 is in a complex environment with strong PAH 
emission. We did not model it in this paper and showed the maps 
and the ratio map for completeness. 

The presented results encourage further exploration of the 
size distribution parameter space. Observational efforts to gain 
information at shorter wavelengths will be important for supply¬ 
ing further constraints. Finally, full 3D modeling based on more 
realistic density structures is required with the option of also fit¬ 
ting thermal emission maps. 
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Appendix A: Determination of the off-core surface 
brightness 


In this appendix, we describe how we have estimated the surface 
brightness near to but off the core and which variation in the de¬ 
termined ratios can be expected when this off-core measurement 
is performed differently. 

The off-core surface brightness needs to subtracted in Eq. ([8]) 


Robs(y, z) 




(A.l) 


since the IRAC measurements are not absolute and contain in¬ 
strumental, as well as back- and foreground, contributions. Since 
4 //(y, y) is needed at any PoSky location of the core, but where 
it cannot be measured, an approximate value needs to be deter¬ 
mined for each point (x,y). Ideally, the region where 4// is mea¬ 
sured is chosen (i) to be near the core (to represent the surface 
brightnesses at core location as closest as possible), (ii) to avoid 
outer core parts (which emission should not be subtracted), (iii) 
to contain no stellar contributions, and (iv) to enable interpola¬ 
tion of surface brightness variations across the core. 
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Fig.A.l. Entire IRAC images containing the core L260 in the 
white frame (see also Fig. 2) for 3.6 yum (top) and 4.5 yum (bot¬ 
tom), respectively. 



In former work, constant Ipff were determined by circular 
averages around the core (e.g.. lNielbock et all2012h or choosing 
the lo cal region near th e core with the lowest surface brightness 
fe.g.. IStutz et al.ll2Q0^ . For cuts through the image, variations 
have also been used, e.g., by linearly int erpolating Ipff from lo¬ 
cations left and right from the core ?e.g.. lSteinacker et ^120101: 
lAndersen et al.ll2Ql^ . 

As is visible from the ex tinction features of cores at 8 yum 
(see, e.g. JPagani et aLll^OlObh . almost all cores have a shape de¬ 
viating from simple spherical symmetry. Judging where outer 
parts of the core gas are located and where we are looking at 
variations in the foreground or background gas is therefore dif¬ 
ficult, especially for cores near the Galactic plane when the FoS 
likely crosses other regions. 

In the following, we use the core F260 with the strongest 
coreshine surface brightness in our sample to illustrate how we 
have performed the off-core measurement and how the results 
depend on the choice of the background. 

In Fig. lA.li the two IRAC images are shown for F260 in the 
bands 3.6yum (top) and 4.5 yum (bottom), respectively. The white 
frame shows the region around the core that is used for the top 


panels in Fig. [51 As is visible from both images, there is a large 
scale horizontal gradient across the background of the core. 

We have chosen three spots near the core to measure G// 
and indicate their location in Panel A of Fig. lA.2l as numbered 
white frames. We also give for comparison the frame 0 in which 
we measured R. In Panel B we show the number distributions 
of pixels as a function of their surface brightness for the entire 
IRAC image. The distributions contain the stellar sources that 
are visible in the high SFB wing of the distribution as they are 
above the mean SFB of each frame: the frame around Core (A), 
the off-frames 1-3, and the R-measurement frame 0. The gradi¬ 
ents visible in the IRAC image results in a mean surface bright¬ 
ness shift from 1 to 3. We therefore interpolate loff from 1 and 3 
which leads to a value close to the mean surface brightness seen 
in Frame 2. 

The situation changes for 4.5 yum. Using the same frames 
as indicated in Panel C, the surface brightness distributions in 
Frames 2 and 3 are almost identical. Nevertheless, the gradient 
between 1 and 3 remains, and we also interpolate loff at 4.5 yum 
from both frames. 

To estimate the uncertainty in the derived range of observ¬ 
able R, we compare the mean R from interpolating between 1 and 
3 and between 2 and 3. Using the mean surface brightnesses of 
the four frames, we get R=2.3 for Framel to 3 off measurement, 
and R=2.41 for Frames 2 to 3. We performed this procedure of 
testing the variation in various frames for all sources discussed 
here. 


Appendix B: Mie calculations 

The absorption and scattering coefficients of large grains can be 
exactly computed for sphe rical particules only using Mie theory 
(iBohren & Huffmanll 19831) . Other grain shapes rely on approx¬ 
imate numerical models, or their validity is restricted to small 
grain sizes. 

In this stu dy the water optical c onstants are extracted from 
the d atabase of Hudgins et al.l (Il993l) . The silicate data are taken 
from iDraine & Feef( 1984. and the amorphous carbon data are 
provided by Zubko et al.l (1 19961) . 

To simulate porous grains, we use an effective medium for¬ 
mulation where the inclusions are made of vacuum. Amorphous 
carbon is also considered as inclusions in the silicate matrix. 

Consider a particulate composite consisting of a matrix 
and including various sizes and shapes made of material other 
than the matrix. Under certain conditions, the composite can 
be homogenized; i.e., the composite can be replaced by a 
homogeneous dielectric medium with the same macroscopic 
electromagnetic response and a certain effe ctive perm i ttivity . 
iFandau & Fifshitzl (119601) and independently iFooveng^ (11964 
proposed an effective medium formulation that take inclusion 
connections for all volume fractionss into account (hereafter 
FFF model). The effective dielectic function is a volume- 
fraction weighted-average of the spherical constituents for the 
composite and is correct to the second order in the differences 
in permittivities, although dipole-dipole interaction is still not 
taken into account. For N constituents, the effective dielectric 
function is 


4"" =< 


1/3 




1/3 


i=l 


(B.l) 


where €{ is the dielectric function of the material /, and ft are the 
volume fractions of arbitrary shape. Here, N is the total num¬ 
ber of inclusions of a given composition. The sum of all volume 
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Fig. A.2. Background choice for L260 


fractions has to be lower than 1. For two constituents, the formu¬ 
lation is extremely simple: 

= (B.2) 

where 6niat is the dielectric function of the matrix. It is clear that 
the LLL model is symmetric with respect to the constituents. 
Other formulations of the effective medium die lectric function 
exist (iMaxwell Gamed[T904 lBruggem^ll935h . each of them 
with their own strengths and weaknesses. We chose to use the 
LLL model for its validity for all volume fractions. 

Water ice is assumed to form a mantle on the top of the 
porous silicate -f amorphous carbon core. Absorption and scat¬ 
tering cross-sections, as well as the phase-functions for ice- 
coated porous spherical grains, are computed using the lat- 
est version of the Mie ro utine for coated spheres provided by 
iToon & Ackermanl(ll98ll) . 
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